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Abstract 

The Persistent Turning Walker Model (PTWM) was introduced by Gautrais et 
al in Mathematical Biology for the modelling of fish motion. It involves a nonlin- 
ear pathwise functional of a non-elliptic hypo-elliptic diffusion. This diffusion solves 
a kinetic Fokker-Planck equation based on an Ornstein-Uhlenbeck Gaussian process. 
The long time "diffusive" behavior of this model was recently studied by Degond 
& Motsch using partial differential equations techniques. This model is however in- 
trinsically probabilistic. In the present paper, wc show how the long time diffusive 
behavior of this model can be essentially recovered and extended by using appropriate 
tools from stochastic analysis. The approach can be adapted to many other kinetic 
"probabilistic" models. 
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1 Introduction 

Different types of models are used in Biology to describe individual displacement. For 
instance, correlated/reinforced random walks are used for the modelling of ant, see e.g. 
[3, 25], and cockroaches, see e.g. [14] and [5] for a review. On the other hand, a lot of 
natural phenomena can be described by kinetic equations and their stochastic counterpart, 
stochastic differential equations. The long time behavior of such models is particularly 
relevant since it captures some "stationary" evolution. Recently, a new model, called the 
Persistent Turning Walker model (PTWM for short), involving a kinetic equation, has 
been introduced to describe the motion of fish [10, 6]. The natural long time behavior of 
this model is "diffusive" and leads asymptotically to a Brownian Motion. 

The diffusive behavior of the PTWM has been obtained in [6] using partial differential 
equations techniques. In the present work, we show how to recover this result by using 
appropriate tools from stochastic processes theory. First, we indicate how the diffusive 
behavior arises naturally as a consequence of the Central Limit Theorem (in fact an 
Invariance Principle). As expected, the asymptotic process is a Brownian Motion in space. 
As a corollary, we recover the result of [6] which appears as a special case where the variance 
of the Brownian Motion can be explicitly computed. We finally extend our main result to 
more general initial conditions. We emphasize that the method used in the present paper 



is not restricted to the original PTWM. In particular, the hypotheses for the convergence 
enables to use more general kinetic models than the original PTWM. 

The present paper is organized as follows: in Section 2, we recall the PTWM and its 
main properties, and we give the main results. Section 3 is dedicated to the proofs. 



2 Main results 

In the PTWM, the motion is described using three variables: position x G M^, velocity 
angle 9 £ M, and curvature k £ M. For some fixed real constant a, the probability 
distribution p{t, x, 9, k) of finding particles at time t in a small neighborhood of (x, 9, k) 
is given by a forward Chapman-Kolmogorov equation 

dtp + T.V^p + Kdep - d^{Kp) - dl^p = (2.1) 

with initial value where 

t{9) = (cos sin 0) = e^^. 
The stochastic transcription of (2.1) is given by the stochastic differential system [t > 0) 

dxt = T{9t) dt 

d9t = Kt dt (2.2) 
dKt = —Ktdt + \/2adBt 

where {Bt)t>o is a standard Brownian Motion on M?. The probability density function 
p{t,x,9, k) of {xt,9t,nt) with a given initial law pQdxd9dn is then solution of (2.1). 
Also, (2.1) is in fact a kinetic Fokker-Planck equation. Note that {Kt)t>o is an Ornstein- 
Uhlenbeck Gaussian process. The formula 



0t = 6Q+ / ds 
Jo 

expresses {9t)t>o as a pathwise linear functional of {Kt)t>o- In particular the process {9t)t>o 
is Gaussian and is thus fully characterized by its initial value, together with its time 
covariance and mean which can be easily computed from the ones of {nt)t>Q conditional 
on ^0 and kq- The process {9t)t>o is not Markov. However, the pair (0j, Kt)t>o is a Markov 
Gaussian diffusion process and can be considered as the solution of a degenerate stochastic 
differential equation, namely the last two equations of the system (2.2). Additionally, the 
process {xt)t>o is an "additive functional" of {9t, i^t)t>o since 

Xt = xq + j T{9s)ds = Xq + j t(^q + j KudiL^ds. (2.3) 

Note that xt is a nonlinear function of {9s)o<s<t due to the nonlinear nature of r, and thus 
{xt)t>o is not Gaussian. The invariant measures of the process {9t,Kt)t>o are multiples 
of the tensor product of the Lebesgue measure on M with the Gaussian law of mean zero 
and variance a^. These measures cannot be normalized into probability laws. Since r is 
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27r-periodic, the process {Ot)t>o acts in the definition of xt only modulo 2tt, and one may 
replace Ohy := M/27rZ. The Markov diffusion process 

{yt)t>o = (^f ^i)t>o 

has state space x M and admits a unique invariant law fi which is the tensor product 
of the uniform law on with the Gaussian law of mean zero and variance , namely 

dn{0,K) = --^=1=151 (^)exp (^-^^dOdK. 

Note that {yt)t>o is ergodic but is not reversible (this is simply due to the fact that the 
dynamics on observables depending only on 6 is not reversible). The famous Birkhoff-von 
Neumann Ergodic Theorem [16, 22, 13, 17, 7] states that for every /i-integrable function 
/ : 5^ X M ^ M and any initial law (i.e. the law of yo), we have, 

pflim (j ff{ys) ds - I fdi^=^= 1. (2.4) 

Beyond this Law of Large Numbers describing for instance the limit of the functional (2.3), 
one can ask for the asymptotic fluctuations, namely the long time behavior as t ^ cxo of 



'^t{\ [ f{ys)ds - [ fdfi 

Jo Js^xR 



(2.5) 



where at is some renormalizing constant such that fif — > 00 as t — > 00. By analogy with the 
Central Limit Theorem (CLT for short) for reversible diffusion processes (see e.g. [9, 17]), 
we may expect, when / is "good enough" and when at = a convergence in distribution 
of (2.5) as t — > 00 to some Gaussian distribution with variance depending on / and on 
the infinitesimal dynamics of {yt)t>o- This is the aim of Theorem 2.6 below, which goes 
actually further by stating a so called Invariance Principle, in other words a CLT for the 
whole process and not only for a fixed single time. 

Theorem 2.6 (Invariance Principle at equilibrium). Assume that yo = {9q,k,q) is dis- 
tributed according to the equilibrium law /.t. Then for any C°° bounded / : S"^ x M ^ M 
with zero fi-mean, the law of the process 

(4)t>o ■= [ fiys)ds, yt/^2] 

\ Jt>o 

converges as e ^ to W-^ (8> where is the law of a Brownian Motion with variance 

Vf = - J gLg dfj. = 2a^ J {df.gl'^ dn 
where L = a^df^ — nd^ + ndg acts on 2-K-periodic functions in 6, and (7 : S*^ x M — > M is 

g{y) = -E ( / f{ys) ds yo = y 
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In other words, for any fixed integer k > 1, any fixed times < ti < ■ ■ ■ < t^, and any 
bounded continuous function F : (M. x S 1 X R)'^ ^ M, we have 

lim E , . . . , 4. )] = E \f{{wI ,Y^), n)) 



where Yi,...,!^ are independent and equally distributed random variables of law ^. and 
where iwi)t>0 ^ Brownian Motion with law , independent ofYi, . . . ,Yk. 

Theorem 2.6 encloses some decorrelation information as e goes to since the limiting 
law is a tensor product (just take for F a tensor product function). Such a convergence 
in law at the level of the processes expresses a so called Invariance Principle. Here the 
Invariant Principle is at equilibrium since yo follows the law fi. The proof of Theorem 2.6 
is probabilistic, and relies on the fact that g solves the Poisson^ equation Lg = f. Note 
that neither the reversible nor the sectorial assumptions of [9] are satisfied here. 

Theorem 2.6 remains valid when / is complex valued (this needs the computation of 
the asymptotic covariance of the real and the imaginary part of /). The hypothesis on 
/ enables to go beyond the original framework of [6]. For instance, we could add the 
following rule in the model: the speed of the fish decreases as the curvature increases. 
Mathematically, this is roughly translated as: 

f{y) = f{e,K) = c{\K\){cose, sine) (2.7) 

where s i— > c{s) is a regular enough decreasing function, see Figure 1 for a simulation. 

The following corollary is obtained from Theorem 2.6 by taking roughly / = r and by 
computing Vr explicitly. In contrast with the function / in Theorem 2.6, the function r 
is complex valued. Also, an additional argument is in fact used in the proof of Corollary 
2.8 to compute the asymptotic covariance of the real and imaginary parts of the additive 
functional based on r (note that this seems to be missing in [6]). 

Corollary 2.8 (Invariance Principle for PTWM at equilibrium). Assume that the initial 
value yo = {9q, kq) is distributed according to the equilibrium fx. Then the law of the process 

Ti9s) ds, yt/e2 (2.9) 

/ t>0 

converges as e ^ to yy^0/U®°° where is the law of a 2- dimensional Brownian Motion 
with covariance matrix ID/2 where 

e-'^'^'-'+'-'Us. 







It can be shown that the constant B which appears in Corollary 2.8 satisfies to 
11 

D = lim — Var(xJ) = lim — Var(xj) where {xj,xf) = xt = I T(9s)ds 

t—*oo t t^oo t Jq 

see e.g. Figure 2. Corollary 2.8 complements a result of Degond & Motsch [6, Theorem 
2.2] which states - in their notations - that the probability density function 

p^t,x,9,K) = ^p 



^It is amusing to remark that "poisson" means "fish" in French. 
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Figure 1: An example of the trajectory t ^ xt = {xl,xf) of the PTWM where the speed 
of the fish decreases with higher curvature (eq. 2.7). Here a = 1 and c{\k\) = 1/(1 + 2|k|). 
The simulation is run during 10 time units, we plot a point each .1 time unit. 
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Figure 2: Convergence of t^^Var(xt) to the constant D. Here a = 1. 



converges as e — > to the probability density 

-^n°(t,x)M(K) 

where M is the Gaussian law with zero mean and variance a^, and solves the equation 

dtn^ - ^ID)A^n° = 

where B is as in Corollary 2.8. Convergence holds in a weak sense in some well chosen 
Banach space, depending on the initial distribution. The meaning of is clear from the 
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stochastic point of view: it is the probabihty density function of the distribution of the 
rescaled process (recall that x is two-dimensional) 

In other words, the main result of [6] captures the asymptotic behavior at fixed time of 
the process (2.9) by stating that for any t, and as e — > 0, the law of this process at time 
t tends to the law of {V^Wt,9, M) where {Wt)t>o, and {6_,M) are independent, {Wt)t>o 
being a standard Brownian Motion, and {9_, M) a random variable following the law /i. 
This result encompasses what is expected by biologists i.e. a "diffusive limiting behavior" . 

Starting from the equilibrium, Corollary 2.8 is on one hand stronger and on the other 
hand weaker than the result of [6] mentioned above. Stronger because it is relative to 
the full law of the process, not only to each marginal law at fixed time t. In particular 
it encompasses covariance type estimates at two different times. Weaker because it is 
concerned with the law and not with the density. For the density at time t we recover 
a weak convergence, while the one obtained in [6] using partial differential equations 
techniques is of strong nature. We should of course go further using what is called "local 
CLTs" , dealing with densities instead of laws, but this will require some estimates which 
are basically the key of the analytic approach used in [6] . 

Our last result concerns the behavior when the initial law is not the invariant law /x. 

Theorem 2.10 (Invariance Principle out of equilibrium). The conclusion of Corollary 2.8 
still holds true when yo = (^c '^o) is distributed according to some law v such that di^sg/dfi 
belongs to L''(^) for some sq > and q > 1, where the law of ys^^. This condition is 

fulfilled for instance if du/dfi belongs to L'^(^) or if v is compactly supported. 

3 Proofs 

The story of CLTs and Invariant Principles for Markov processes is quite intricate and it 
is out of reach to give a short account of the literature. The reader may find however a 
survey on some aspects in e.g. [9, 15, 27, 17]. Instead we shall exhibit some peculiarities 
of our model that make the long time study an (almost) original problem. The under- 
lying diffusion process {Ot,Kt)t>o with state space M? is degenerate in the sense that its 
infinitesimal generator 

L = a^dl2 - Kd^ + Kde (3.1) 

is not elliptic. Fortunately, the operator dt + L is Hormander hypo-elliptic since the 
"diffusion" vector field X = (0, a^) and the Lie bracket [X, Y] = XY - YX of X with the 
"drift" vector field Y = (k, — k) generate the full tangent space at each {9,k) G M^. The 
drift vector field Y is always required, so that the generator is "fully degenerate". This 
degeneracy of L has two annoying consequences: 

1. any invariant measure i/ of L is not symmetric, i.e. J fLgdv ^ J gLf dv for some 
nice functions / and g in L^(z^), for instance only depending on 6. 

2. the carre du champ of L given here by Tf = ^L{f'^) — fLf = 2q:^|c?k/P is degenerate, 
so that one cannot expect to use any usual functional inequality such as the Poincare 
inequality (i.e. spectral gap, see [1, 23]) in order to study the long time behavior. 
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This situation is typical for kinetic models. In the more general framework of homoge- 
nization, a slightly more general version of this model has been studied in [11], see also 
the trilogy [19, 20, 21] for similar results from which one can derive the result in [6]. The 
main ingredient of the proof of Theorem 2.6 is the control of the "rate of convergence" 
to equilibrium in the Ergodic Theorem (2.4), for the process {6^, nt)t>o- We begin with a 
simple lemma which expresses the propagation of chaos as e goes to 0. 

Lemma 3.2 (Propagation of chaos). Assume that yo = {9jq,ko) is distributed according to 
the equilibrium law ii. Then the law of the process {y^)t>Q = {yt/e'^\yQ converges as e — > 
to //®°°. In other words, for any fixed integer k > 1, any fixed times < ti < • • • < tk, 
and any bounded continuous function F : {S^ X M)'= R, we have 

limE[F{y!^,. . . ,yfj] = E[FiY,, . . . ,n)] 

where Yi, . . . ,1^^ are independent and equally distributed random variables of law /i. 

Proof of Lemma 3.2. Let us denote by L the operator (3.1) acting this time on 27r-periodic 
functions in 9, i.e. on functions S*^ x R — > M. This operator L generates a non-negative 
contraction semi-group {Pt)t>o = (e*'^)t>o in L^(/i) with the stochastic representation 
Ptf{y) = IE[/(ys)|?/o = y] for all bounded /. We denote by L* the adjoint of L in L^(/-i) 
generating the adjoint semi-group P^ , i.e. 

L* = a^9^ — nd^ — ndg 

acting again on the same functions. The function H{y) = H{6, k) = 1 + satisfies 

L*H = -2H + 2(a2 + 1)<-H + 2{a^ + 1) l|,|<^/2^ (3.3) 

so is a Lyapunov function in the sense of [2, Def. 1.1]. Since C = S^x{\k\ < V2^^^^TT} 
is compact and the process {yt)t>o is regular enough, C is a "petite set" in the terminology 
[2, Def. 1.1] of Meyn & Tweedie. Accordingly we may apply [2, Th. 2.1] and conclude 
that there exists a constant K2 > such that for all bounded / satisfying J fd^ = 0, 

||Pi/|lL2(^) < K2 \\f\\^e-\ (3.4) 

We shall give a proof of the Lemma for k = 2, the general case k > 2 being heavier but 
entirely similar. We set s = ti < t2 = t. It is enough to show that for every bounded 
continuous functions F,G:S^ x M ^ R, we have the convergence 

hm E[F(yf)G(yf )] = E[F(y)]E[G(y)] 

where y is a random variable of law /i. Since yo follows the law /i, we can safely assume 
that the functions F and G have zero //-mean, and reduce the problem to show that 

E[F{yt)G{yt)] = J P,,AFPit-s)/e^G) dfi = J FP(^t_,y,2G dfi — ^ 0. 
Now since /i is a probability measure, we have L^(^) C l^^{fJ^) and thus 



- WF{t-s)/e^G\\^ < ||FP(j_s)/e2G||2 < \\F\\oo\\P{t-s)/e'^G 



7 



The desired result follows then from the 



L°° bound (3.4) since 



{t-s)/e'- 



■ 2 G 



L2(m) 



<if2||G|Le-(*-^)/^' — >0. 



□ 



Proof of Theorem 2. 6. The strategy is the usual one based on Ito's formula, Poisson equa- 
tion, and a martingale CLT. However, each step involves some peculiar properties of the 
stochastic process. For convenience we split the proof into small parts with titles. 

Poisson equation 

Let L, L*, and {Pt)t>o be as in the proof of Lemma 3.2. Since / is bounded and satisfies 
J fd^ = (i.e. / has zero /i-mean), the bound (3.4) ensures that 



It follows that g belongs to the L^(/_f)-domain of L and satisfies to the Poisson equation: 

Lg = f ■ml.\^l). 

Since /i has an everywhere positive density with respect to the Lebesgue measure on xM, 
we immediately deduce that g belongs to the set of Schwartz distributions V and satisfies 
Lg = / in this set. Since L is hypo-elliptic (it satisfies the Hormander brackets condition) 
and / is C°°, it follows that g belongs to C°° . Hence we have solved the Poisson equation 
Lg = f in a strong sense. Remark that since g G ]L^(^) and / is bounded, we get 




Furthermore, the formula Ptg — g = /q Pgf ds ensures that 



Jim - (^t - 5) = / strongly in {p) . 



E^b/] < 00. 



Ito's formula 



Since g is smooth, we may use Ito's formula to get 




Lg{ys)ds almost surely 



which can be rewritten thanks to the Poisson equation Lg = f as 




almost surely. 



(3.5) 



This last equation (3.5) reduces the CLT for the process 
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to showing that {£{g{yt/e^) ~ 9{yo)))t>o goes to zero as e — > and to a CLT for the process 



aeV2 J d^,g{ys)dBs 



For such, we shall use the initial conditions and a martingale argument respectively. 
Initial condition 

Since the law /i of yo is stationary, Markov's inequality gives for any constant K > 0, 

n\9{yt/e^)\ > K/e) = n\9{yo)\ > K/e) < ^ 0. 

It follows that any n-uple of increments 

- 9{yto), ■ ■ ■ ,9{ytr,) - 9{yt„-r)) 

converges to in probability as e ^ 0. Thanks to (3.5), this reduces the CLT for 



e / f{ys)ds 
. ^0 



to the CLT for 



Martingale argument 

It turns out that ((Mf )j>o)e>o is a family of local martingales. These local martingales 
are actually martingales whose brackets (increasing processes) 



Jo 

converge almost surely to 

2aHEi,[\d^g\'^] = tVf as e ^ 

thanks to the Ergodic Theorem (2.4). According to the CLT for L^-martingales due to 
Rebolledo, see for example [12] for an elementary proof, it follows that the family {Mf)t>o 
converges weakly (for the Skorohod topology) to Vj {Bl)t>o where {Bl)t>o is a standard 
Brownian Motion. Consequently, we obtain the desired CLT for the process 



e / f{ys)ds 

Jo . 



t>0 



Namely, its increments are converging in distribution as e — > to the law of a Brownian 
Motion with variance Vr- It remains to obtain the desired CLT for the process (-zf)f>o- 
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Coupling with propagation of chaos and asymptotic independence 

By the result above and Lemma 3.2, the CLT for (^f )j>q will follow if we show that 

'0 



f{ys)ds\ and {yt/e^)t>0 = {It/e-^^ l^t/e^) 
■^0 /t>o 



t>0 



are independent processes as e — > 0. It suffices to establish the independence as e — > for 
an arbitrary /c-uple of times = to < < ■ ■ ■ < ^fc = ^- To this end, let us introduce a 
bounded continuous function h and the smooth bounded functions 



where 1 < j < A; for given real numbers 6i , . . . , 6^ . Let us define 

k / 



A' = E„ 



Kytle-2)WhAe 

j=i V 



fiVs) ds 



Introduce = {t/e^) — {t/y/e) and = t/y/e. For e small enough, > {tk-ie^), so that 



using the Markov property at time we get 



A' = E,, 



WhAe 

3=1 V "-1 



nys)dsjE^ 

The conditional expectation in the right hand side is equal to 



hklel f{ys)ds]E 



h{ys,)hk{ £ I f{ys)ds 



yo = yu 



and the second term can be replaced by 



E 



Kys,] 



yo = yt. 



up to an error less than 



going to as e — > 0. It thus remains to study 
k-l / 



E„ 



'J 



f{ys)ds\hk{e I f{ys)ds\E 



KysA 



yo = yu 



Conditioning by yt^ , this can be written in the form 

H{e,y)PsMy)^^{dy) 



with a bounded H, so that using the convergence of the semi-group, we may again replace 
Ps^h by / hdjj, up to an error term going to 0. It remains to apply the previously obtained 
CLT in order to conclude to the convergence and asymptotic independence. □ 
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Remark 3.6 (More general models). The proof of Theorem 2.6 immediately extends to 
more general cases. The main point is to prove that g solves the Poisson equation in L^(^). 
In particular it is enough to have an estimate of the form 

\\Ptf\\i.^,)<(^{t) 
for every t > with a function a satisfying 



a{s) ds < oo. 

According to [2], a sufficient condition for this to hold is to find a smooth increasing 
positive concave function ip such that the function a defined by 

1 1 

ait) = 1 where G^{u) = / — 7-^ ds 

satisfies the integrability condition above, and a Lyapunov function H > 1 such that 

J Hdn< 00 and L*H < -ip{H) + 0{lc) 

for some compact subset C. In particular we may replace the Ornstein-Uhlenbeck dynam- 
ics for K by a more general Kolmogorov diffusion dynamics of the form 

dKt = -VV{Kt)dt + V2dBt. 

The invariant measure of {Kt,Ot)t>o is then e~^^'^^dKd9. We refer for instance to [8, 2] for 
the construction of Lyapunov functions in this very general situation. For example, in one 
dimension, one can take V'{x) = |x|^ for large and < p < 1. Choosing H{y) = \k\'^ 
for large k. furnishes a polynomial decay of any order by taking q as large as necessary. 
Actually, in this last situation, the decay rate is sub-exponential, see for example [8, 2]. 

Remark 3.7 (Asymptotic covariance). It is worth noticing that if the asymptotic variance 



{AV)f = }im jE, 



t \ 2- 

f{ys)ds 







exists, then Vf = {AV)f. Similarly we may consider complex valued functions / and 
replace the asymptotic variance by the asymptotic covariance matrix which takes into 
account the variances and the covariance of the real and imaginary parts of /. 

Proof of Corollary 2. 8. We may now apply the previous theorem and the previous remark 
to the 2-dimensional smooth and /i-centered function r. The only thing we have to do is 
to compute the asymptotic covariance matrix. To this end, first remark that elementary 
Gaussian computations furnishes the following explicit expressions 

Kt = HQ + V2a [ e^^^dBs, (3.8) 
Jo 

et = eo + {l- e-*) Ko + V2a [ (1 - e^"*) dB, . (3.9) 

Jo 
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Since xj = cos 6 gds and x^ = sin 6 gds, Markov's property and stationarity yield 

ft 



[xj 



E. 



E„ 



(xl sin 9s + x1 cos ( 



ds 



(3.10) 



(cos 6u sin 9s + sin 9u cos 6s) du ds 



Jo 



[ [ E^[sm{9u + es)]duds 
Jo Jo 

sin [ 290 + 2(1 - e"-^)Ko + V2a 



JO 



(l_e^-(«-«))dS; 



ds 



where (i?j)t>o is a Brownian Motion independent of (kqj^o)- Since kq and ^^'^ ^^^o 
independent (recall that ;U is a product law), we may first integrate with respect to 9q 
(fixing the other variables), i.e. we have to calculate E^(sin(2^Q + C)) which is equal to 
since the law fi of 9q is uniform on [0, 27r[. Hence the /i-covariance of {xl,x^) is equal to 
(since this is a Gaussian process, both variables are actually independent), and similar 
computations show that the asymptotic covariance matrix is thus D/2 where 



□ 



Proof of Theorem 2.10. We assume now that yo ~ ^ instead of yo ~ ^- We may mimic the 
proof of Theorem 2.6, provided we are able to control Ey{g^{ys)). Indeed the invariance 
principle for the local martingales (Mf )t>o is still true for the finite-dimensional conver- 
gence in law, according for instance to [13, Th. 3.6 p. 470]. The Ergodic Theorem ensures 
the convergence of the brackets. The first remark is that these controls are required only 
for s > So > where sq is fixed but arbitrary. Indeed since r is bounded, the quantity 



TiOs) ds 







goes to almost surely, so that we only have to deal with J^^^ 



ISO 



so that we may replace 

by So in all the previous derivation. Thanks to the Markov property we thus have to 
control Ejy^^ (5^(2/5)) for all s > 0, where i^so denote the law of yso- This remark allows 
us to reduce the problem to initial laws which are absolutely continuous with respect to 
/i. Indeed thanks to the hypo-ellipticity of ^ -|- L we know that for each sq > 0, Vsq is 
absolutely continuous with respect to fi. Hence we have to control terms of the form 



E„ 



di', 



so 



d/j, 



{yo)9^{ys) 



The next remark is that [2, Theorem 2.1] immediately extends to the framework for 
2 < p < CO, i.e. there exists a constant Kp such that for all bounded / satisfying J fdfi = 0, 



(3.11) 



Since the function / is bounded and satisfies J fdfj, = 0, the previous bound ensures that 
g belongs to LP(/Lt), for all p < 00. In particular, as soon as dvso/dfi belongs to E''(//) for 
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some I < q, g{ys) belongs to L^(Pi,) for all s > sq. Additionally, these bounds allow to 
show without much efforts that the "propagation of chaos" of Lemma 3.2 still holds when 
the initial law is such a u. To conclude we thus only have to find sufficient condition for 
dvso/dpi to belong to one IL'^(^) {q > 1) for some sq > 0. Of course, a first situation is 
when this holds for sq = 0. But there are many other situations. 

Indeed recall that for non-degenerate Gaussian laws r/i and r/2 the density dr]2/dr]i 
is bounded as soon as the covariance matrix of r/i dominates (in the sense of quadratic 
forms) the one of r]2 at infinity, i.e. the associated quadratic forms satisfy qi{y) > q2{y) 
for \y\ large enough. According to (3.8) and (3.9) the joint law of (K.t,6t) starting from a 
point (k, 0) is a 2-dimensional Gaussian law with mean 

mt = {e~^K,e + (1 - e"*)K) 

and covariance matrix Dt = c? At with 

/I - (1 - e-^f \ 

^* " V(l - e-*)2 2t - 3 + 4e-* - e'^* ) ' 

Note that if the asymptotic covariance of (kj, dt) is not 0, the asymptotic correlation van- 
ishes, explaining the asymptotic "decorrelation" of both variables. It is then not difficult 
to see that \i v = by is a Dirac mass, then dug/d^i is bounded for every s > 0. Indeed for 
t small enough. At is close to the null matrix, hence dominated by the identity matrix. 
It follows that dvt/dr] is bounded, where ?7 is a Gaussian variable with covariance matrix 
The result follows by taking the projection of 6 onto the unit circle. A simple 
continuity argument shows that the same hold if is a compactly supported measure. □ 

Remark 3.12. Once obtained such a convergence theorem we may ask about explicit 
bounds (concentration bounds) in the spirit of [4] (some bounds are actually contained in 
this paper). One can also ask about Edgeworth expansions etc. However, our aim was just 
to give an idea of the stochastic methods than can be used for models like the PTWM. 

Remark 3.13. The most difficult point was to obtain L^(/i) estimates for d^Q- Specialists 
of hypo-elliptic partial differential equations will certainly obtain the result by proving 
quantitative versions of Hormander's estimates (holding on compact subsets U): 

\\d.g\\^<C{U){\\g\\^+\\Lg\\^). 

We end up the present paper by mentioning an interesting and probably difficult 
direction of research, which consists in the study of the long time behavior of interact- 
ing copies of PTWM-like processes, leading to some kind of kinetic hypo-elliptic mean- 
field/exchangeable Mac Kean-Vlasov equations (see for example [24, 18] and references 
therein for some aspects). At the Biological level, the study of the collective behavior at 
equilibrium of a group of interacting individuals is particularly interesting, see for instance 
[26]. 
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